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Abstract: The geodesic problem in Wasserstein spaces with a metric 
perturbed by a conformal factor is considered, and necessary optimahty 
conditions are estabihshed in a case where this conformal factor favours 
the spreading of the probability measure along the curve. These conditions 
have the form of a system of PDEs of the kind of the compressible Euler 
equations. Moreover, self-similar solutions to this system are discussed. 

1 Introduction 

Let us consider a closed convex set Q C M*^ and the set V{Q) of probability 
measures in il. Given p G (l,oo), we denote by Wp{Q) the subspace of 
measures with finite p-th moments, i.e. 

Wp{n) := |;u G V{n) : J \x\Pdfi < cx)| . 

We endow Wp{i^) with the canonical Wasserstein distance Wp{fi, v) of order 
p (see P, for the basic facts about Wp). 

It is well known that Wp{^) is a length space, and that (constant speed) 
geodesies of Wp(0) are in one to one correspondence with optimal transport 
plans, via McCann's linear interpolation procedure (see for instance Propo- 
sition 7.2.2 of ^). Here we consider, instead, the case when the Wasserstein 
metric is perturbed by a conformal factor L[^): by minimizing 

\\^it)\^l'Wt)dt (1.1) 





among all curves /i connecting fiQ = fi to fii = , one obtains a new squared 
distance depending on p and L, and we are interested in computing the 
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geodesies relative to this distance. In is the rate of change of 

Wp also called metric derivative, along the curve ^, see 1)2. 

This problem has been introduced in where the main goal was 
to choose a factor L favouring atomic measures in order to give a time- 
dependent approach to some branched transport problems which may be 
applied to the study of river networks, pipe systems, blood vessels, tree struc- 
tures. ... In fact, by setting L{fi) = (for < r < 1) if = ^ - aiSx^ 
and L = +oo on measures which are not purely atomic, there is a strong 
link between this variational problem and those which were first presented 
in and [HI (the latter uses in fact a time-dependent approach, but by 
means of measures on the space of paths instead of paths on the space of 
measures). This choice of L is in fact a local functional on measures which, 
among probability measures, favours the most concentrated ones. In [2], as 
a natural counterpart, the case of local functionals L which prefer spread 
measures is considered as well and the two problems sound somehow specu- 
lar. The aim of the present paper is in fact to consider this second problem 
and to find out optimality conditions in the form of PDEs. 

In particular, we study in detail the case when L(^) is the 7-th power 
of the L"? norm of the density of // with respect to Lebesgue measure C^, 
with q > 1 and 7 > given, and L{fi) = +00 if is a singular measure. 
So, geodesies with respect to the new metric tend to spread the density as 
much as possible. Denoting by ut the density of fit, we find that a necessary 
optimality condition for geodesies is (for p = 2, see ()2.5() for general p) 

{K{t)vu) + K{t)V ■ {v O vu) + H{t)Vu'i = 0, (1.2) 

where Vt is the tangent velocity field of /Uj, linked to ut via the continuity 
equation + V • (vtUt) = 0. Here H{t) < and K{t) > are suit- 
able functions depending only on the metric derivative of fit and on L{fit). 
As Brenier pointed to us, this equation is very similar to the compressible 
Euler equation, but with a negative pressure field p = H{t)u'^; a similar 
equation, with H constant and q = 3, recently appeared also in in the 
one-dimensional case. In fact the main difference appears in the relationship 
between the L part and the speed part: here it is multiplicative, while in |7j 
it is additive, as we will explain in a while. 

The appearence of the Euler equation as an optimality condition is not 
very surprising, taking into account the approach developed, in the incom- 
pressible case, by Brenier (first in a purely Lagrangian framework in [1], 
and then in a mixed Eulerian-Lagrangian one in |^, [Sj). In this connec- 
tion, we mention that our derivation of the optimality condition differs from 
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[1], [S], where duality is used to perform first variations, and uses instead a 
perturbation argument directly at the level of the primal problem. 

Due to the non-convex nature of this problem, we don't know of any 
sufficient minimality condition for the geodesies. In this connection, one 
may notice that, in the case 7 = g/2 and p = 2, we have 



inf (5 / u''dx + l / Ivl'^udx = ILiuC^) ( / \v\ 
s>o Jn S Jn \Jn 



"^udx 



1/2 



and the minimal norm of v is strictly linked to the metric derivative. 

This suggests a connection between the "multiplicative" model studied here 
and in [21, and the "additive" model 



mm 



{[ [ u'' + \v\'^udxdt : ^u + V-{vu)=0 
[Jo Jn dt 



subject to Dirichlet conditions at t = and t = 1. This additive model, 
in the case g = 3, is exactly the one studied in [7] (in this connection, see 
also [S]). Notice that this problem is convex in the pair (u,vu). It turns 
out, indeed, that the (necessary and sufficient, by the convex nature of the 
problem) optimality conditions for the additive model are very similar to 
(|1.2|1 . the only difference being that H and K do not depend on time. 

In the last part of the paper we compute and characterize particular 
self-similar or solutions of 1)1.2(1 . 

Acknowledgements. We warmly thank Y. Brenier for many useful com- 
ments on an earlier version of this paper and G. Buttazzo for the interest 
towards this work, which in particular lead us to study the self similar so- 
lutions of Section 3. 



2 Optimality Conditions for Weighted 
Wasserstein Geodesies 

2.1 A new velocity vector field 

Definition 1. If we are given a Lipschitz curve ^ : [0,1] — > Wp{il.), we 
define velocity field of the curve any vector field : [0, 1] x ^ M"^ such 
that for a.e. t G [0,1] the vector field Vt = v{t,-) belongs to [LP(/Ut)]'^ and 
the continuity equation 

d , . 

-/it + V • (yiit) = 
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is satisfied in the sense of distributions: this means that for all (p G C^{^) 
and any ti < t2 £ [0, 1] it holds 

(pdfit^ - f (pdfit^ = [ ds [ Vcf) ■ Vgdju^s, 



ti Jn 



or, equivalently, in differential form: 
d 



dt 



/(t)dfit= / ^(t)-Vtd^t for a.e. te[0, 1]. 
Jo. 



We say that v is the tangent field to the cm've /it if, for a.e. Vt has minimal 
[U'{^t)]'^ norm for any t among all the velocity fields. 

It is now well known (see for instance Theorem 8.3.1 and Proposi- 
tion 8.4.5 in [T]) that for any Lipschitz or absolutely continuous curve /ij 
with values in Wp{Vt) there exists a unique tangent field and moreover it is 
characterized by 

The right hand side, in the equality above is the rate of change of Wp along 
the curve /it, also called metric derivative of /if. 

We want now to investigate how velocity fields change if we modify the 
curve /if. 

Theorem 2.1. Let a Lipschitz function /if : [0, 1] — > Wp(r2) and a smooth 
function T : [0, 1] x ^ 6e given, such that for any t the function 
Tf := T{t, •) is a difjeomorphism. Let us consider the new curve fi[ given by 
= (Tt)f^nt- If Vt is a velocity field for fit, then the vector field v' defined 

by 

" dT\ ' 

VTf • 'i;f + — /if 



v't- l^'t = {Tt)i 

is a velocity field for 
Proof. We have 



dt 



^ M+h - 4'dHt= 4>o Tt+hdfit+h - I (poTtdfit 
Q JQ JQ Jn 

= {(f>o Tt+h - o Tf ) dut+h + / (1)0 Tt d{nt+h - l^t) 
JQ Jn 

In {It ^^'^^ ° ■ '^^) dfit+h+J^ ds J^{V(t>)oTfVTfVs dfi^ 
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where in the last equahty we have used the fact that Vt is a velocity field 
for n, with test function (j) o Tf. It is now convenient to divide by h, rewrite 
and pass to the limit as /i — > 0: 

+ / (V(/)) o Tt ■ VTt -Vidiit + T I ds I ViJt ■ (vsd^s - Vtdnt) , (2.2) 
Jn Jt J(i 

where ipt = o Tf. In the first term on the right hand side the measures 
^lt+h weakly converge to ^f, since t i-^ /i^ is Lipschitz continuous, while 
the integrand uniformly converges as a function of the space variable x to 
(V(?!>) °Tt ■ ^ as /i — 0. Hence we get convergence of the integral. If we 
prove that the last term tends to zero at least for a.e. t G [0, 1] we get the 
thesis, since then we would have 

lim ^^^t+h ~ 
h 

= [{S4>) °^f^ + (V'^) ° ■ VTt ■ vt^ dut = J^^'t'- vM, 

and this is nothing but the differential version of the continuity equation 
for v' and n' (it remains to prove v'^ G L'P{ii[) but this is straightforward 
since Tt is a diffeomorphism and this allows to write down the densities and 
estimate them). To prove that the last term vanishes at the limit we see 
that, for fixed il) G Lip(J7) the function 



gAs) := [ Vi^- Vgdjis = T" / i^'^l^s 
Jn ds Jq 



is L°° since fit is a Lipschitz curve in Wp{n) and hence almost any s G [0, 1] 
is a Lebesgue point. This allows to fix a negligible set C [0, 1] such that 
any point t G [0, 1] \ A'^ is a Lebesgue point for all the functions g^^_^ for 
ti G Q. We fix now t G [0, 1] \ A'' and try to prove that the last integral in 
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(tO) tends to zero. For ti G Q it holds 



1 f 

-J ds J Vtpt ■ {vsdfis - vtdfit 



< 



t+h 



ds 



+ 



V(V't -ipti) ■ Vtdfit 



+ 



^ r-t+h r- 

-J ds J Vtpti ■ {VsdfJLs - Vtdfit 



t+h 



ds / VV'ti • {vgdus - vtd^t 



In the last sum the second term tends to zero by the fact that i is a Lebesgue 
point for g^^_^ and the first term may be made as small as we want by choosing 
ti close to t, since ipg = (j) ° Ts and both cp and T are regular. □ 



2.2 Derivation of the optimality conditions 

We consider the minimization problem presented in [^j, i.e. finding a curve 
of measures in Wp{fl) of minimal length according to a metric which, roughly 
speaking is the Wasserstein (infinitesimal) metric multiplied by a conformal 
factor. Precisely, if we define for q > 1 the functional 

1 +00 otherwise, 
we want to minimize ^ 

/ Lg{fit)\fi'mdt, 

Jo 

where is the metric derivative of the curve fj, and the minimization 

occurs among all the Wasserstein-Lipschitz curves t ^ fit with given initial 
and final points, i.e. fiQ and fii are given probability measures in Wp(0). 
We will always consider only the non trivial case fi,Q ^ fii. If we define 
V{fi,t) = \vt\^dfit, where v is the tangent field to the curve fit, we know 
that \fi'\{t) = V{fi,ty/P. We may generalize the functional we want to 
minimize by considering 

difi) := [' L,{fitrV{fi,tfdt 
Jo 

which reduces to the case studied in [2j if a = 1 and (3 = 1 /p. Notice that in 
this case the functional does not change under reparametrization of curves. 
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while if (3 > 1/p the minimization selects a particular parametrization. For 
/3 < 1/p the existence of a minimum is not ensured. Anyway we do not 
deal here with existence results (see but we only look for necessary 
optimality conditions. We will consider variations of ^ of the form 

= (Tf with T%t,x)=x + e^{t,x), = id + e^{t,-), 

for arbitrary regular functions ^ G C^([0, 1] x r2;R'^). In the end optimality 
conditions will be expressed through a system of PDEs: we will obtain the 
result after collecting some lemmas. What we want to do now is exploit- 
ing the fact that for a minimizing curve ^ the following quantity must be 
minimal for e = 0: 

^{,jii) = l^j\,{trv,{tfdty 

provided we define -^^(i) = Lq{nl) and Ve{t) = V{^^,t). Since it is not 
completely easy to deal with the term it) , we will replace it by (t) , with 
Ve{t) given by 

VS)= [ \{vnt\'df,t. 
Jn 

Here the vector field is the one we get by Theorem 12 . 1 1 when the map T 
is given by and the initial field Vt is the tangent field to fXf In this way 
we have Vs{t) > V^(i) (since vf is not necessarily of minimal norm) but 
^^(O) = 14(0). Thus we may switch to considering 1^(0 instead of Vs{t), 
getting 

We will compute the derivative of 5^(//f) with respect to e and get the con- 
ditions we are looking for. 

Lemma 2.2. If jj is a curve given by fit = and such that 5^(/i) < +oo, 
then for almost any t G [0, 1] it holds 

In particular, if we compute the derivative at e = 0, we have 
j-FS)\e=^ = {l-q) jjyOuldC''. 
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Moreover, for e sufficiently small ( depending on T, hut not ont) the follow- 
ing inequality holds: 

Proof. We look at the integrand function in the definition of F^: to do this it 
is necessary to look at the density of the measure juf . Thanks to the change 
of variables formula, this density can be easily seen to be given by 



6 /T^e 



where J stands for the Jacobian (this formula is a consequence of being 
a diffeomorphism at least for small e). Thus, after changing variables, we 
have 

The derivative of the integral is given by 

Ut 



(l-g)(jr, 



where {JTl)' stands for the derivative w.r.t. e of JT^. This quantity may 
be easily estimated by Cuf, since 1 — a < JTf < 1 + a and {JTf)' < B 
for suitable constants a and B. Since for almost any t the function ut must 
belong to L'^ (because the functional we are minimizing is finite) we can 
apply the dominated convergence theorem and get the thesis. To obtain the 
derivative at £ = it is sufficient to notice that {JTfy\s=o = V • ^, which is 
well-known. The same estimate we used to get dominated convergence may 
be used to get the last inequality. □ 

In the next lemma we consider the term V^. 

Lemma 2.3. If fx is a curve such that ^{fJ-) < +oo, then for almost any 
t G [0, 1] it holds 



(2.3) 

In particular, if we compute the derivative at e = 0, we have 
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Moreover, for e sufficiently small ( depending on T, hut not on t) the follow- 
ing inequality holds: 

^V,it)<CiVifi,t) + l). 
de 

Proof. If we compute the densities of /if and the expression of the new 
velocity field and we change variable in the integral by y = Tf{x), as we did 
in the previous lemma, we get 



Veit) 



VTf • vt + 



dt 



dfit 



(2.4) 



When we differentiate the integrand we get exactly the integrand in (|2.3|) . 
and we need only to show that this expression is uniformly dominated, at 
least for small £ and almost every t to get the result. By boundedness of the 
derivatives of T"^ it is not difficult to see that the norm of the first vector in 
the scalar product in the integrand may be estimated by 



VTf • Vt + 



while for the second it holds 



dt 



p-i 



< {C\vt\ + C) 



< C\vt\+C 



p-i 



for a suitable constant C. Hence, since Vt G [LP{fit)Y for almost every t the 
integrability is proved and the differentiation under the integral sign can be 
performed. □ 

To conclude, we must put together the two previous results in order to 
compute the derivative of the integral in t. 

Theorem 2.4. If ^ is a curve with 5^(/i) < +c>o and V{fi, t) > Vq > for 
almost every t, then it holds 



^S(.')l„o 



a{l-q) Tf^-V^ / {V-Ou'ldC^dt 
Jo Ja 

+ PP ^ .vt + ^)-Vt dfit dt, 



where F{t) = Lq(fj,t) CLnd V{t) has the usual meaning. 
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Proof. By the definition of we see that the pointwise derivative of 

the integrand is given by aFs{t)"-^^Ve{tf + /3Fe(t)"K(i)'^"^^- By the 
regularity of the term Fi;{t) may be estimated both from above and below 
by F{t), up to multiplicative constants. As far as V^{t) is concerned, the 
argument is a little bit more tricky. Indeed we must write V^{t) according 
to (|2.4|) . then estimate 

VTF ■vt + 



A'\vt\ - B< 



< A^\vt\ + B, 



dt 

for e small enough, where the constants are as close to 1 as we want and 
the constant B is as small as we want (this comes from VT^*" = id -\- 0(e) 
and arVat = 0(e)), and get 

A-y^ -B <V^ < + B. 

The assumption y > Vq > allows us to infer from these inequalities 
that also may be estimated both from above and below by V up to 
multiplicative constants. Finally, by the estimates in Lemmas 12.21 and 12.31 
we bound the whole pointwise derivative by CF'^V^ since we had 

^ < OF; ^ < C{V + 1) < 0(1 + Uv, 
de de Vq 

where the last inequality too comes from V >Vo. Since L'^V^ is integrable 
on [0, 1] , we may differentiate under the integral sign and get 

^5(^^)|,=o = ^' (^aF{tr-'^U=oV{tf + pF{trV{tf-'^U=o^ dt. 

The result follows when we replace the derivatives in e by the explicit ex- 
pressions we computed in Lemmas 12 . 21 and 12 . 31 □ 

Remark 1. If /? = 1/p and is a minimizer, it is always possible to get the 
lower bound V > Vq hy reparametrizing in time, for instance by choosing 
the constant speed parametrization. 

Corollary 2.5. If fi minimizes ^ with given boundary conditions /Uq and 
fii, then its density u and its tangent field v satisfy 

a{l-q) C F{t)''~^V{tf [ {V ■ Ou't dC^ dt 
Jo Jn 

+ p(3 [ F{trV{tf-^ [ ut\vtr\V^-vt + ^)-vtdC''dt = 0, 
Jo Jn <Ji 

for any vector field ^ G O;?°(]0, l[xJ]; ] 
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Proof. It is sufficient to notice that when we create the modified curve fj,^ 
starting form the vector field ^ we do not change the initial and final points 
of the curve, so that the minimality implies that the derivative of ^{^J''^) at 
e = vanishes. □ 

2.3 A system of PDEs 

The following theorem follows directly from the previous section. 

Theorem 2.6. Let /xq, fJ-i G Wp{^}) and let fi be a curve minimizing ^ on 
r(/io,/Ui), with a finite minimum value. Then, denoting byu{t,-) the density 
of nt and by v{t,-) the tangent field to the curve fi, {u,v) provide a weak 
(distributional) solution of the system 

'H{t)Vu'i + K{t)V ■ {u\v\P''^v ®v) + f^ {K{t)u\v\P-'^v) =0 in 17, 
jlU + V ■ {vu) =0 inVL 
uv ■ n = on d^l 

lim u(t, AC^ = liQ] limn(t, = ^1, 

(2.5) 

where H{t) = a(l - q)F{t)'^-^V{tf and K{t) = p(3F{tYV{tf-^ . 
Given (/Uq, ^i), existence of minimizers is ensured whenever q <l + \/d or, 
for general q, under the assumption that fiQ = uqC^ and fii = uiC^ with 
uq, ui E L'^{Q) (see f^j), hence under these conditions existence of solutions 
to this system is ensured. 

It is interesting to rewrite the equations, make some formal simplification 
and look at some particular cases. 

First we expand all the terms in the first equation of System 1)2. 5() . ob- 
taining 

H{t)Vu'' + K{t) {u\v\P~^v ■ Vv + v\v\P-'^V ■ {uv) + u{v- V\vY'-'^) v) 

+ K(t) ( v\v\P-^^u + u4- iv\v\P-^)] + —K{t)u\v\P~^v = 0. (2.6) 
\ dt dt ^ J dt 

Notice that this is always a vector equation, i.e. a system itself, consisting of 
d equations with d+ \ unknown functions (the components of v and the den- 
sity u). This system is then completed by the continuity equation. As usual, 
by V ■ Vv we mean the vector whose i— th component is Y^j{vjdvi/dxj). 

A formal simplification in (|2.6|) may be done: in fact there is a term 
{K{t)v\v\P~'^){du/dt + V • (uv)) that might be removed by using the conti- 
nuity equation. This is actually possible only under extra regularity assump- 
tions on K and v (it consists of testing the continuity equation against the 
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product K{t)v\v\P ^ which is not in general or regular enough). Anyway, 
after this formal simplification, H2.6|) becomes 

H{t)Vu'^ + K{t) {u\v\P-'^v ■Vv + u{v- V\v\P^'^) v) 

+ K{t)u^ {v\v\P~^) + ^K{t)u\v\P-\ = 0. (2.7) 

Notice that in the case /3 = 1/p we can reparametrize in time the solution 
and there are several possible parametrization choices that present some 
advantages. For instance, we could choose a parametrization so that K{t) 
is constant, to get rid of the final derivative in time. This choice implies 

and this, in the case of a bounded \^\ < +00, is sufficient to have the lower 
bound V > Vq, since in this case F is bounded from below by a positive 
constant. 

Another important fact to be noticed is that in 1)2. 7(1 there is a common 
u factor. It is still formal, but in this way we should get, on {u > 0}, 

H{t)u'''^Vu + K{t) {\v\P''^v -Vv + lv- V\v\P-'^) v) 
+K{t)^ {v\v\P-^) + ^K{t)\v\P~\ = 0. 

(JjL (JjL 

Remark 2. One might wonder whether the solutions u are automatically 
positive a.e. in Q. for t g]0, 1[. This could be suggested by the fact that 
in the minimization problem spreadness of the density is favoured. In next 
session we will see with explict examples that this is not necessarily the case. 

We finish this overview of simplifications of the system by looking at the 
simplest case, i.e. p = q = 2, a = l, j3 = 1/2, in the parametrization regime 
where K is constant. In this case we get 

' -2V{t)^'^Vu + K{v-Vv + f^v) =0 in {n > 0}, 

^u + V ■ (vu) = in 17 

uv ■ n = on dO, (2-^) 
limn(t, •)£'^ = /io; limu(t, •)£'^ = /ii. 
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Under no constraint on the parametrization we have, instead, 



2V{t)^/^Vu + K{t) {v-Vv + j-^v) + = in {-u > 0}, 

(2.9) 



■^u + V ■ (vu) = in 



dt 

uv ■ n = on dQ 

tio tTi 

3 Self-similar solutions 

3.1 Homothetic solutions with fixed center 

In this section we look for particular solutions of the System (|'2.5() which 
are self-similar in the sense that, for any t, the measure fit is the image 
under an homothety of a fixed measure. For simplicity we will consider only 
the case of System ()2.9p . i.e. with p = q = 2, and we assume that G il. 
The regularity of the candidate solutions we will propose will be enough to 
ensure that we can use this simplified system, instead of System (|2.5p . To 
start this analysis it is necessary to establish the following Lemma. 

Lemma 3.1. If fi is a curve in W2{^) of the form fit = (^i?(t))ttA fof « 
certain regular function R : [0, 1] — >]0, 1] (where Tr{x) = Rx is the multipli- 
cation by a factor R, hence an homothety), then its tangent field is given by 
vt{x) = xR!{t)/R{t). 

Proof. It is not difficult to prove that the field we defined solves the conti- 
nuity equation and hence is a velocity field. Indeed, if </> G C^(il), it holds 

A [ ^dfit = ^ I 4){R{t)x) dii{x) = I V 4){R{t)x) ■ R' {t)x dfi{x) 
dt Jq dt Jq Jq 



f V(P{R{t)x) ■ ^^R{t)xdfi{x) = [ V4>-vtdnt. 
Jn Ryt) Jn 



It remains to prove that v is actually the tangent velocity field, i.e. that its 
norm is minimal for a.e. t. This is achieved if we are able to prove that 
ll'^^t = l^l'(i) for a.e. t S [0, 1]. To do this, let us fix two times t < t + h 
and see that the map T[x) = xR{t + h)/R{t) is a transport between /i^ and 
fit+h- Since it is the gradient of the convex function x i— > x'^R{t + h)/2R{t), 
it is actually the optimal transport according to the quadratic cost. Hence 



X diJ,t{x). 



13 



Since this last quantity is exactly the norm of Vt in 
V is the tangent field to the curve /i. 



(//j), this proves that 



□ 



Remark 3. In the case p ^ 2 the same result is true, but one has to use the 
characterization of tangent velocity fields in terms of closure of gradients of 
smooth maps, see Proposition 8.4.5 of 

A first result we prove is the following: 

Theorem 3.2. If {u,v) is a self-similar solution of the system 1)2. 5|) with u 
Lipschitz continuous, then necessarily u is of the form 

u{t,x) = {At — V for suitable coefficients At, Bt > 0. 

Proof. We look at the equation (|2.8|) with p = q = 2, which is valid on 
{u > 0}, and we freeze time, i.e. we look at the resulting space equation 
for fixed t. We use the fact that v is of the form Vt{x) = ctx, which implies 
that any term v, v ■ Vv and dv/dt are of the same form. This easily implies 
that also Vu is of the same form. Hence, at time t, on {u > 0}, it holds 
u{x) = At — Btx"^, where a priori Bt could also be negative. Anyway we can 
prove that Bt cannot be negative. In this case in fact, if were a convex 
unbounded domain, then u could not be the density of a probability measure. 
On the other hand one can easily see that on bounded convex domains 
self-similar solutions must vanish on dQ, otherwise we should get a jump of 
the density at the boundary of {u > 0} when rescaling, but u was supposed 
to be Lipschitz (except in the case that the solution is constant in time). 
This implies that also in the case of a bounded O the coefficient Bt must 
be positive. For the same continuity reason we get that the region {u > 0} 
must agree with the region Q n {At — Btx^ > 0} in order to have continuity 
of u, and this proves the formula. □ 

Remark 4. A similar result could be obtained for generic Wasserstein spaces 
with exponent p > 1, getting that any self-similar solution should be of the 



form u{t, x) = {At - Bt\x\P) V 0. 



Theorem 3.3. If fi is a probability measure on 0, with density 



u{x) = A[{R'^ 



xp) VO] 



then for any regular and monotone function R : [0, 1] [0, 1] the curve 
Ht = (TR(t))[j/x is a solution to the System 1)2. 5() together with its tangent 
field V. 
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Proof. It is sufficient to check the first vector equation in the system ()2.9|) . 
First we compute the correct constant A: we must have 



'd + 2' 

and hence A = R^'^~'^{d+ 2)/{2uJd)- This allows us to compute the term 

2{d + 2) 



1 = A - r^)duJdr'^-^dr = AR 

Jo 

tience A 

F{t): 

F = A^ (R^- r'^fdoJdr'^'^dr = R', . 
Jo {d + ^pd 

Then we compute V by recalling that Vt{x) = xR'{t)/R{t). It holds 
V = (^^"j Aj\\R^ - r')du;dr'-'dr = ^ 

We must also compute dv/dt and v ■ Vv: 

dv R"R-{R'f ^ (R'\r ^ (R'V 

We compute now 

m = F{t)v{t)-'i' = R''\R'\-'^=^, 

^Jd{d + A)ujd 

d-l r>-dr r>'\-2 r>r'\ 2(d + 2) 



K'{t) = sign(R'){-dR-'^-^ - R-'^{R')-^R 



y/d{d + A)iOd 



If we can c = sign(R')^M±^ it holds K = cR-'^(R')-^ and K' = 

c{-dR-'^-^ - R-'^iR'y^R"), but also -2V^/'^Vu{x) = cdR'R-'^-^x. In- 
serting everything in the equation we must check that 

dR'xR-'^-^ + R-\r')~^x^ - (dR-'^-^ + R-'^{R'y^R")x% = 0. 

it R 

The proof is achieved as this last equation is (miracolously enough) always 
satisfied. □ 

Remark 5. By a similar proof we can show that, for p 7^ 2, if /2 has a density 
of the form u{x) = A[{RP — \x\^) V 0], then ^ gives raise to a self-similar 
solution. 

Remark 6. This kind of self-similar solutions can join two different prob- 
ability measures which are homothetic, and in particular arrive up to the 
Dirac mass 5q. Anyway it is not in general possible to link a measure to 5q 
by a curve with finite energy: in j^, conditions to ensure this possibility are 
provided, but in general they are not satisfied in the case q = 2. 
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3.2 Moving self-similar solutions 

We have characterized all the self-similar solutions which link two homo- 
thetic probability measures. It is however interesting to look also at the 
moving self-similar solutions, i.e. at solutions obtained by homotheties and 
translations together. 

In this case we consider a reference measure fl and we look for solutions 
of the form {T^)^fl, where T*(x) = R{t)x + x{t). It is not difficult to replace 
Lemma l3 . 1 1 with the following: 

Lemma 3.4. If fi is a curve of the form fit = (2^*)tlA; then its tangent field 
is given by 

vt{x) = ^{x-x{t))+x'{t). 

Proof. The result may be proved very similarly to Lemma l3.ll it is sufficient 
to check the continuity equation 



4>{R{t)x + x{t)) dnix) = / V(j){R{t)x + x{t))-{R'{t)x + x'{t))dfi{x) 



d_ 

R' it] r 

V<i){R{t)x + x{t)) ■ -^{R{t)x + x{t)) dfi{x) = J^VcP- vt dfit, 
and then to check the optimality of the norm by the fact that the map 
X ^ ^^^^'^^^^\ x - x{t)) + x{t + h) 

transports nt on fit+h and is optimal, and that 

1 f (R{t + h) , , , V , , . 

y j^^^^ [x - x{t)) + x{t + h)-xj dfit{x) 

converges to 

(R'it^ \^ 
-^{X-X{t))+X'{t)j dflt{x) = \\Vt\\L2(^^,,). □ 

For computational simplicity we consider moving self-similar solutions 
only under a special reparametrization. 

Theorem 3.5. If Jl is a probability measure on 0, with density 

u{x) =A[{R^- VO] 
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and x(0), x(l) G 0, are assigned, a curve = (^*)tl/^; parametrised so that 
K = FV^^/"^ is constant, is a moving self-similar solution (solving System 
H2.8|) together with its own tangent field) if and only if the vector x moves 
on the straight line segment from x{Q) to x{\) with constant speed and R is 
a strictly concave function of t. This means 

x" = 0; R?'^{d(yR')'^ + {d + 4)(x')^) is constant and R strictly concave. 

Proof. We only need to check under which conditions the first equation is 
satisfied. We re- write in this case the quantity considered in Theorem 13.31 
first we compute 

n(x) = A[(/?2_|^_^|2)vO]; A = ^^^, ^^^^^ = "T^^^ " 

F = R-'^:^^±^- v = -^(R!f + {x'f- 

We have used the fact that ut is symmetric around x{t) and hence there is 
no mixed term [x — x{t)) ■ x'{t) in computing V{t). Then we go on with 
dv/dt and v ■ Vv: 

dv , ^R"R-{R'f ,R' „ ^ fR'\r 
^ fR'V, , R' _, dv ^ , _R" _„ 

V ■ \/V = [ ^ ] [X — X) + —X + ' vt> = [x — X)— + X . 

\ R J R ot R 

Then we look at the the condition to have K'{t) = 0, which is equivalent to 
F~'^V being constant, and thus R'^'^{d{R')'^ + ((i + 4)(x')^) must be constant. 
Assuming K to be constant we try to satisfy the equation, and we write it 
in the following form that we can reach after multiplying by V^^'^: 

fdv 1 \ 
-2VVu + F [ — + -V - Vv] =0. 
\ot 2 J 

This equation becomes 



d + V ' ' ' ) ujs,R'^^'^ {d^A)u:d 'R 

To satisfy this equation it is necessary and sufficient that the two parts, the 
one involving x — x and the other with x" both vanish. After simplifying we 
get 

R-\d{R'f + {d + A){x'f) + ^ = 0; x" = 0. 

R 
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Hence we must have x{t) = (1 — t)x{0) +tx{l) and x'{t) = e = x{l) — x{0). 

Now we recall that R'^'^{d{R')'^ + {d+A){x')^) was assumed to be constant and 
so d{R')'^ + {d + A){x'f = CR-^'^. Hence we get R" = -CR'^'^-^. Thus, 
M is a moving self-similar solutions if and only if the following conditions 
simultaneously hold: 

'd{R')^ + {d + 4)e2 = CR'^'^ for a certain C, 
< R" = -Ci?-2<^-i for the same C, 

^x{t) = x{0)+te. 

By differentiating the first equation we get 2dR'R" = -2dC R-'^'^'^ R' and 
hence the second is automatically satisfied, provided we can ensure that 
R' ^ a.e. This means that R being strict concave is sufficient (it is not 
possible to have more than a time where R' vanishes) , but it is also necessary 
from the second equation. The result is then proved. □ 
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